Verifying black hole orbits with gravitational spectroscopy 
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Gravitational waves from test masses bound to geodesic orbits of rotating black holes are simu- 
lated, using Teukolsky's black hole perturbation formalism, for about ten thousand generic orbital 
configurations. Each binary radiates power exclusively in modes with frequencies that are integer- 
linear-combinations of the orbit's three fundamental frequencies. General spectral properties are 
found with a survey of orbits about a black hole taken to be rotating at 80% of the maximal spin. 
The orbital eccentricity is varied from 0.1 to 0.9. Inclination ranges from 20° to 160°, and comes 
to within 20° of polar. Semilatus rectum is varied from 1.2 to 3 times the value at the innermost 
stable circular orbits. The following general spectral properties are found: (i) 99% of the radiated 
power is typically carried by a few hundred modes, and at most by about a thousand modes, (ii) 
the dominant frequencies can be grouped into a small number of families defined by fixing two of 
the three integer frequency multipliers, and (iii) the specifics of these trends can be qualitatively 
inferred from the geometry of the orbit under consideration. Detections using triperiodic analytic 
templates modeled on these general properties would constitute a verification of radiation from an 
adiabatic sequence of black hole orbits and would recover the evolution of the fundamental orbital 
frequencies. In an analogy with ordinary spectroscopy, this would compare to observing the Bohr 
model's atomic hydrogen spectrum without being able to rule out alternative atomic theories or 
nuclei. The suitability of such a detection technique is demonstrated using snapshots computed at 
12-hour intervals throughout the last three years before merger of a kludged inspiral. The system 
chosen is typical of those thought to occur in galactic nuclei, and to be observable with space-based 
gravitational wave detectors like LISA. Because of circularization, the number of excited modes 
decreases as the binary evolves. A hypothetical detection algorithm that tracks mode families dom- 
inating the first 12 hours of the inspiral would capture 98% of the total power over the remaining 
three years. 

PACS numbers: 04.70.-s, 97.60. Lf 



I. INTRODUCTION 

The birth of physics in its modern form can arguably 
be placed at the first successful efforts to monitor the mo- 
tion of the planets as they orbit the sun, and to model 
their planar elliptic orbits. Since then, gravitation has 
transitioned from the best to the least tested among the 
fundamental laws of physics. Gravitational spectroscopy, 
monitoring power spectra of the outputs from gravita- 
tional wave detectors, can verify the existence of two- 
body systems with a dramatically different character. 
Athough characterized by mass ratios similar to those 
for planetary systems, these systems are so distorted by 
the strong gravity near black hole event horizons that 
their orbits look more like balls of twine than planar el- 
lipses. These orbits differ significantly from those with 
relativistic precessions which have already been observed 
and which may soon be measured in our galactic center 
Given the relative paucity of gravitational experi- 
ments, and the enigmatic status of gravitation compared 
to the other fundamental interactions, opportunities to 
observe these systems are highly valuable. Though un- 
likely candidates for Brahc-Kcpler-like catalysts of an- 
other revolution in physics, such observations would help 
to elevate tests of gravity to a level that better compares 
with those for the other fundamental physical interac- 
tions. 

General relativity predicts that black holes will radiate 



gravitational waves characterized by discrete frequency 
spectra in at least two broad classes of realistically ob- 
servable scenarios: black hole ringdowns, and during ex- 
treme mass ratio inspirals (EMRIs). In both, the ra- 
diation is a consequence of a single black hole having 
been slightly perturbed. Observations of this radiation 
are therefore ideally suited for studying phenomena gov- 
erned completely by the physics of isolated black holes. 

In the case of black hole ringdown, the source of the 
black hole's perturbation is an arbitrary transient event, 
and the radiation is emitted as the hole settles back into a 
quiescent state. Ringdown radiation consists of a super- 
position of exponentially damped monochromatic waves, 
with frequencies and damping timescales that are de- 
termined by a pair of continuous parameters, the black 
hole's mass and spin. Potentially observable ringdown 
events include binary systems that merge to form sin- 
gle black holes, and supernovae that result in black hole 
formation. Ringdown phenomena have been studied in 
detail elsewhere [1, 0, 0; Q , and are not the focus of this 
work. 

This work instead focuses on EMRIs — black hole bi- 
nary systems in which the mass of the black hole is much 
greater than the mass of the companion. The planned 
space-based interferometer LISA is anticipated to observe 
anywhere from tens to thousands of EMRIs produced by 
the capture of compact stellar mass black holes (and the 
occasional white dwarf) by the megamassive black holes 
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found in galactic nuclei (masses ranging from 10^ Mq to 
10^ Mq) with sufhcient sensitivity to determine the mass, 
spin, and quadrupole moment of the larger black hole to 
within a fraction of a percent 0, H, M, [HI • In the 
recent mock LISA data challenge, a proof of principal for 
a variety of detection algorithms using simulated signals 
and noise, three groups using independent detection al- 
gorithms recovered EMRI masses and spins to within a 
few to a few tenths of a percent [I^, [H, ■ 

Advanced ground-based detectors (with a lower fre- 
quency cutoff of about 10 Hz [l^) might also observe sim- 
ilar capture events by intermediate black holes in globular 
clusters [3 ■ Observation of these intermediate mass ra- 
tio inspirals (IMRIs) are more speculative, might occur 
with an estimated upper limit of a few to a few tens per 
year [l^ , and due to their less extreme mass ratios, may 
require theoretical waveform models that are more so- 
phisticated than those based on black hole perturbation 
theory alone. See Ref. [l^ for a review of EMRI science 
with both ground-based and space-based detectors. 

Unlike radiation from black-hole ringdown, for EMRIs 
the source of the black hole perturbation is persistent. 
On timcscalcs that are short compared to the observable 
lifetime of the binary, the companion behaves as a test 
mass orbiting on a geodesic of the background spacetime. 
Such orbits are in general characterized by three funda- 
mental frequencies [Tgl. [20|. and the corresponding radi- 
ation is tri periodic, or rather, a superposition of modes 
that oscillate exclusively at frequencies that are integer- 
linear combinations of the orbit's three fundamental fre- 
quencies 
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mkn 

OJmkn = mUJ^ + kiOg + ULOr , 



(1.1a) 
(1.1b) 



where I am using the convention that all unconstrained 
summation indices range from —oo to oo. Here w^, we, 
and Wr, are the fundamental orbital frequencies asso- 
ciated with motion in the (Boyer-Lindquist) coordinate 
shown as a subscript, and are the two independent 
components of the metric perturbation measured by a 
distant observer, and hmkn are complex amplitudes that 
depend on the observer's position, and on the parameters 
of the binary (masses, spin, and orbit geometry). 

The radiation described by Eq. (jl.ip is only a snap- 
shot of a complete EMRI waveform. While work toward 
precision waveforms for non-test-mass motion through- 
out the lifetime of EMRIs remains an active field (see 
the discussion of Capra waveforms in Ref. [2^1. for short 
enough times these waveform snapshots are as exact as 
any envisioned. For example, for the case of a mass /i 
on a circular orbit of initial radius r about a nonrotat- 
ing hole with mass M fi, the normalized overlap of 
the radiation from geodesic motion alone and that from 
the true inspiraling motion will be greater than 95% for 



times shorter than [22 



/ M \ / r \ii/4 fM/n\ 
tdephasc « (1 day) ^-^^ j [—) j 



whereas the orbital period is 



'-orbit 



minutes) 
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(1.2) 
(1.3) 



Understanding the radiation from simple test mass mo- 
tion may prove a sufficient basis for useful observations, 
however crude. This is the motivating principal behind 
the work described here. 

This paper has two main results. The first is an ob- 
servation of trends in a survey of numerically simulated 
EMRI-snapshots, for thousands of different generic or- 
bital configurations, using the code described in Refs. [l^ 
and [l^l. The observed trends are summarized as follows: 
(i) When the modes are sorted in order of decreasing 
power, their power decreases fast enough to be primarily 
confined to a relatively small number of modes compared 
to the number of modes computed when using the algo- 
rithm introduced in Refs. [2^ and [2^. (ii) The dominant 
modes can be grouped into a small number of families de- 
fined by fixing two of the integer frequency multipliers, 
(iii) The specifics of these trends (for example, the precise 
falloff of power as a function of mode index, or the iden- 
tity of the dominant mode families) can be qualitatively 
inferred from the geometry of the orbit under considera- 
tion, as illustrated by the following two examples. 

Simple orbits that are very nearly circular are expected 
to be typical of IMRIs that could be found with ground- 
based detectors. Mandcl et al. estimate that the most 
likely IMRI formation mechanism should result in an ec- 
centricity e < 10^'' by the time the orbital frequency is 
brought into the observable band above 10 Hz for ad- 
vanced LIGO [iBl, and that even for the less likely for- 
mations mechanisms that can yield higher eccentricities, 
90% of the systems should have e < 0.1 at 10 Hz [13. 
For these orbits power falls off as a power law with mode 
index, and 99% of the radiated power is confined to a few 
to ~ 10 modes with frequencies (jl.lbp in the following 
families: 

■muj(f, dominant power at m = 2 

uj = i ztmujij, + ujr dominant power at m = 2 . 
ztmu!^ + oje dominant power at m = 1 

(1.4) 

where the upper (lower) sign of ± refers to prograde (ret- 
rograde) orbits, and mode-power decreases exponentially 
as TO is varied away from its value for the dominant mode. 
Note that for power spectra, there is no observational 
consequence for a change in a frequency's overall sign. 
For any given spectrum, I always define mode indices 
(to, k, n) in such a way that the frequency is positive: 

^mkn > 0. 

More eccentric orbits are expected to be typical of the 
EMRIs observed by space-based detectors. EMRIs that 
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fall into LISA'S frequency band are thought to be born 
with such large mitial eccentricities (10^^ < 1— e < 10"'^) 
that, although radiation circularizes them, even at the 
time of merger as many as half of the systems should have 
residual eccentricities of about e > 0.2 [l3|. Increasing 
orbital eccentricity dramatically slows the rate of power 
falloff and ultimately results in a spectrum dominated, 
albeit much less so, by the following mode families: 

{±2w0 + nojr 
±UJ^ + LUg + nujr ■ (1.5) 
zt2uj0 + nuij. 

where again ± is + for prograde orbits and — for ret- 
rograde orbits, and mode power decreases as n is varied 
away from its value nmax at which the spectra are peaked. 
The value of rimax can be very crudely predicted, with an 
accuracy on the order of 10%, from eccentricity e accord- 
ing to the formula 

rimax « exp(l/2) (1 - e)-^/' , (1.6) 

that was fitted [IJ] to the peaks in the spectra derived 
by Peters and Mathews, who in 1963 used an analytical 
treatment of Newtonian orbits with the quadrupole for- 
mula for gravitational radiation [2^. For these eccentric 
orbits, the rate of power falloff as a function of mode 
index is such that as many as ^ 10^ to 10^ modes are 
needed to capture 99% of the power radiated. 

The second main result of the paper is a proposal for 
combining the analytical waveform model (jl.ip with the 
observed spectral trends so as to create a means for ver- 
ifying fundamental aspects of black hole physics. Preci- 
sion measurements that will test the fundamental physics 
of black holes have long been among the primary motiva- 
tions for gravitational wave observatories [2^, [23| . How- 
ever, whether or not there are any feasible means for im- 
plementing one of the more grand measurements, mea- 
suring the multipolar decomposition of spacetime near 
black hole candidates in such a way as to identify gen- 
eral relativity as the only valid theory, has been a subject 
of debate [1^, ■ In its most pure form, that measure- 
ment will require more than the current understanding 
of what exactly is predicted by general relativity and 
its alternatives. I propose a related, but less ambitious 
measurement that will constitute a minimal verification 
of one of general relativity's fundamental predictions for 
black hole perturbation: that black holes perturbed by 
bound test masses radiate according to the above wave- 
form model (|l.ip for times that are long compared to 
the spectrum's fundamental periods 27r/uj^^e,r but short 
compared to the inspiral time. The measurement would 
also recover the evolutionary sequence of those three fun- 
damental frequencies. This is a minimal verification in 
that it only provides a means for confirming general rel- 
ativity's prediction under the assumption that the per- 
turbed spacetime geometry is that of a rotating black 
hole. Observations that provided such a verification may 
ultimately be subjected to the more grand tests based 



on generalizations [s^l of Ryan's theorem (sTj. In the 
absence of a greater theoretical understanding of alter- 
natives to general relativity and black holes , however, the 
verification alone would not disqualify alternatives for ei- 
ther the background spacetime or the theory of gravity. 

The remaining sections of this paper are outlined as 
follows. In Sec. |TT1 I review the relevant equations that 
define generic black hole orbits and that describe the ra- 
diation they produce. In Sec. Illli the radiative spec- 
tra for three sample configurations are studied and gen- 
eral trends are discussed. In Sec. IIVI I simulate spectra 
from a grid of orbital configurations and discuss how the 
general trends from the previous section vary across the 
grid. Section |V] describes how the spectral trends evolve 
throughout the lifetime of an approximated inspiral. In 
Sec. I VII I outline a practical means for extracting sig- 
nals characterized by these general spectral trends from 
data collected by gravitational wave observatories. Sec- 
tion |VlT] summarizes the paper's main points. 



II. EMRI-SNAPSHOT SPECTRA 

Radiation from generic configurations of test masses 
bound to black holes has been studied in previous work 
[H, I'^ this section, I review definitions and equa- 
tions derived there which will be needed throughout the 
remaining sections. 



A. Orbits 

The physical system described by an EMRI waveform 
snapshot is a nonspinning^ test mass /i bound to a black 
hole with mass M and a spin per unit mass of magnitude 
< a < M. The test mass' orbit is a bound geodesic of 
the Kerr spacetime determined by M and a. In Boyer- 
Lindquist coordinates {t,r,9,(f>), the orbit as a function 
of proper time (r) satisfies the four first order geodesic 
equations derived by Carter [H, H^l . When bound solu- 
tions are parametrized as functions of Mino time A, with 
dr = (r^ -|- cos^ 9)dX, two of their coordinates are pe- 



For the EMRIs that could be seen with LISA, the mass ratio 
renders the spin of the smaller object negligible (see Appendix C 
of Ref . [Toll ) . The same may not be true for the IMRIs that could 
be seen with advanced ground-based detectors, but the effect 
would likely be competing with more significant complications 
due to the less extreme mass ratios for IMRIs. 
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riodic [1^ 

<(A) ^r\ + J2 tfe„e-*(^'^«+"^'-)^ , (2.1a) 

kn 

r(A) = ^r„e-™^'-\ (2.1b) 

n 

0(A) ^^efce-*"^''^ (2.1c) 



(/)(A) = T^A + ^ 0fc„ 



(2.1d) 



where the coefficients in front of the exponentials are 
constants (with values that cause these seemingly com- 
plex sums to be real). The quantity F relates the Mino- 
frequencies T^fi^ to coordinate-time frequencies 



(2.2) 



that appear in the radiation observed by distant ob- 
servers (jl.ip . The three spatial Boyer-Lindquist coordi- 
nates of the orbit are not periodic functions of coordinate 
time i, however it follows from the formalism in Ref. [2lj 
that they have simple biperiodic forms 



r{t) = E 

kn 



'^kn^ 



kn^ 



-i(kLO0-\-nLJr)t 



-i{ku}Q-\-niJ.r)t 



(j){t) ^cj^t + y^^fc^e" 

kn 



(2.3a) 
(2.3b) 
(2.3c) 



where the expansion coefficients are again constants. A 
derivation of the coefficients ffc„, and (jjkn in terms 
of the coefficients in the Mino-timc series (|2.ip is given 
in Appendix [XI 

The orbital frequencies uJr^e.tp are uniquely deter- 
mined by specifying the three constants associated with 
Killing fields, energy E, axial angular momentum L, and 
Carter's constant Q 



E ~ —f^iu , 
L = fiMu'l' , 

+ COt^ e + a2(^2 _ ^2-^ ^Qg2 ^ 



(2.4) 
(2.5) 

(2.6) 



where u" = dx'^''/dT is the orbit's four-velocity. They 
can also be determined by specifying the orbit's coordi- 
nate boundaries between two radial turning points and 
between two angular turning points that arc symmetric 
about the equatorial plane at 7r/2 



^min — ^ ^ ^max 5 



(2.7) 
(2.8) 



or by specifying three geometric constants generalized 
from Newtonian orbits: eccentricity e, semilatus rectum 



p, and inclination t 



sgn(a;0)6'n 



P 



M 



M 



1 + e 
P 



1 



(2.9) 
(2.10) 
(2.11) 



where sgn(a;0) is 1 for prograde orbits and —1 for ret- 
rograde orbits. See Appendix A of Refs. [2^ or [2^ for 
explicit formulae relating the geometric orbital constants, 
the formal Killing constants E, L, and Q, the frequencies 
^4>.e.r^ and the Mino frequencies T and T^^e.r- Each of 
the following sets of parameters are uniquely determined 
by fixing the values for any one of them 



(T^,Te,TO , 

(^miii; ^maxi ^min) 

(e,p,i) . 



(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 



Following the terminology of the Guelph group [35|, [Sy, 
[srj . each of these triples is a complete set of principal 
orbital elements. 

After fixing the principal orbital elements, the orbit is 
not completely determined until one specifies an initial 
position, or some equivalent set of parameters, called po- 
sitional orbital elements [3^, [3^, [331 ■ Here I will use the 
following orbital elements 



(At, A0, A,., Xg) 



(2.17) 



These are defined such that, after specifying them and 
the principal orbital elements, any bound black hole orbit 
can be uniquely expressed as follows 



t(A) 



r(A) 

e{\) 

m 



V{\-Xt) + Y,iisin[kTe{X-\e)] 

k=l 

t;sin[7iT^(A- A,.)] , (2.18a) 

Ef„cos[nTr(A- A^)] , (2.18b) 

^feCos[A:Te(A- Ae)] , (2.18c) 



E 



n=0 
oo 

E 

k=0 



T^{X~X^) + Y,^l sin[fcTe(A 

fe=i 



A«)] 



<f>j^ sin[?iT,.(A — Xr)] 



(2.18d) 



were the hatted coefficients depend only on the principal 
orbital elements and arc given by integrals over the A- 
derivativcs of the coordinates as given by Eqs. (2.27) and 
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(2.28) in Rcf. [2J|^. The first two positional elements, At 
and can take on any value. The second two are defined 
on < Ar < 27r/Tr and < Ae < 27r/Te, respectively as 
the smallest positive values of A at which the coordinate 
shown in their subscript reaches the smaller of its two 
turning points. 



B. Radiation 

Black holes that are forever perturbed by bound test 
masses produce a gravitational wave field whose two or- 
thogonal linearly polarized components, /i+ and hx , can 
be expressed as a single complex function made up of 
a series of modes that oscillate at frequencies that are 
integer-linear combinations of the orbit's fundamental 
frequencies [2l|. For an observer located at (i, r, 0, 0), 
that function is 



h+ — ihx = 



mkn 



up to corrections of order /i^/M^. 
amplitudes arc given by 



(2.19) 



The complex mode 



^rakn 



ruj 



1=2 



id) (2.20) 



i{9) and Zimkn are quantities that are found 
by a ppl ying Teukolsky's black hole perturbation formal- 
ism |38| . as described in Sec. Ill of Ref. [IJ. The real 
fimction Simkn{(^), that has here been normalized over 
the unit sphere 



at radial infinity, the complex cocfRcicnts i?°,j*fc„ deter- 
mine the amount of outgoing radiation (toward infinity) , 
while the Rp^kn determine the amount of in-going radi- 
ation (toward the event horizon). The complex numbers 
Zimkn appearing in the mode amplitudes for the grav- 
itational wave (|2.20p are taken from the solution that 
obeys a boundary condition with zero ingoing radiation 
at infinity 



R, 



imki 



(2.23) 



as well as only in-going radiation at the event horizon. 
Details on how to calculate these solutions, and on how 
the horizon-boundary condition enters, can be found in 
Ref. [13. 

The time-averaged rates of change for the principal or- 
bital elements can be given in terms of the quantities 
described above'^. Each rate of change can be expressed 
as two fluxes: an outgoing flux at radial infinity and an 
ingoing flux at the event horizon. Horizon fluxes are 
usually two or three orders of magnitude smaller than 
the fluxes at infinity (see Rcf. [i^ for an analytic treat- 
ment of circular nonspinning binaries, or Tables IV and 
VII from Ref. [l^l for numerical examples of generic bi- 
naries). Though they may prove observable indirectly 
through thcir influence on the orbital evolution [30| . I 
will not discuss them here apart from mentioning that 
their calculation is nearly identical to that of fluxes at 
infinity. 

The time-averaged power radiated to infinity at fre- 
quency LOmkn, as measured by distant observers, is given 
by m 



27r / dO [Shnkn{0)f sine ^ I 



(2.21) 



satisfies Teukolsky's angular equation, a homogeneous 
ordinary differential equation. The complex numbers 
Zirnkn are constant coefficients that come from the lim- 
iting behavior of the physical solutions to Teukolsky's 
radial equation, an inhomogeneous ordinary differential 
equation. At infinity, that equation's general solutions 
have the following form (see Teukolsky's original papers 
[s^ . or for a particularly good textbook- level discussion. 
Sec. 4.8.1 of [3^) 



DOUt 

^hnkn 



R 



1 



Imkn , 



Rlmkn (r OO) 

(2.22) 

The solution of Teukolsky's master equation fully de- 
scribes the system's radiation, to first order in the pertur- 
bation, and is a sum over products of the separated func- 
tions Rirnknir), Si,nkn{9), 6™"^, and e"'"-'"'*. Therefore, 



dE 



rnkn 



Attlu 



mkn 



■El 

1=2 



(2.24) 



and the total time-averaged energy flux at infinity is then 
just a sum of the power radiated at all possible frequen- 
cies 



dE 
'dt 



E 

mkn 



dE 



(2.25) 



mkn 



The remainder of this paper studies the distribution of 
power among the various terms in this sum and how that 
distribution is affected by the configuration of the orbit. 
Understanding this distribution is akin to understanding 
the evolution of the other princple orbital elements, since 
they arc all somewhat simply related. The time-averaged 
flux of the other two formal constants of motion can be 



^ The time averages can be simply defined in terms of Mino-time 
^ Because of a difference in notation, replace Ax there with x used integrals of length 27r/Tr and 2-K/lg. See Sec. 3.8 and 9.1 of 

here, for x = t,<h. Ref. 1331 for details. 
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written as [H, [H, [HI, [H 




where the angled brackets represent a time average. From 
these expressions, one can compute the timc-averaged 
evolution for any set of principal orbital elements due to 
radiation at infinity. 

This all but concludes the review of quantities needed 
here to describe EMRI snapshots. Before moving on how- 
ever, an important property about how the various or- 
bital parameters influence the radiation should be ad- 
dressed. It has been shown that the positional orbital 
elements have a somewhat simple influence on the ra- 
diation, by comparison to the influence of the principal 
orbital elements. Their only influence on the quantities 
discussed above is as a phase factor in the complex num- 
bers Zimkn ■ That factor can be written as [H, [43l | 

Zimkn{\t,\,p,K,>^e) = e'^'-'-ZimUO, 0,0,0) , (2.28) 

where Xmkn is given by 

Xrnkn = TTlT ^{X^ - At) + kT g{\g - Aj) + nT,.(A,. - Aj) . 

(2.29) 

The radiative fluxes for the principal orbital elements 
(|2.26p are therefore independent of the positional orbital 
elements, since they depend only on \ Zimkn\- So the posi- 
tional orbital elements will not be relevant to the discus- 
sion of power spectra and their evolution in the remaining 
sections. However, the mode amplitudes of the waveform 
(j2.20p are functions of Zimkn, and not just their moduli. 
Therefore, knowing the evolution of the principle orbital 
elements alone is insufficient for evolving from one snap- 
shot to the next in an optimal coherent matched filtering 
detection algorithm. To evolve waveforms coherently, one 
needs a prescription for changing both the principal and 
positional orbital elements. This issue will be revisited 
when discussing detection algorithms in Sec. IVII 

Before discussing the results of the simulations, it is 
perhaps useful to remind readers who are more familiar 
with other simulations of radiating black hole binary sys- 
tems that for these snapshot spectra, the source of the 
radiation is ever-present. There is no initial data from 
which imperfections could produce junk radiation which 
dies out over time. Of course no source of radiation could 
really persist forever like this, but that is why these spec- 
tra are "snapshots" of EMRI spectra. As described in the 
introduction, the snapshot spectra match the true spec- 
tra from EMRIs only for sufficiently short observation 
times, shorter than the dephasing time (|1.2p . 



III. SAMPLE SPECTRA 

In this section I describe EMRI snapshot spectra from 
three sample systems and identify properties that are use- 
ful for understanding spectra from systems with arbitrary 
orbit geometries. All of the spectra shown in this paper 
were simulated u sing the numerical code that was first 
described in Ref. [2J|. When simulating the spectra de- 
scribed in this paper, all adjustable parameters of that 
code were set to the same values used when computing 
the catalog of orbits introduced in Sec. V of that paper, 
with one exception. The one exceptional code parameter 
is the requested fractional accuracy £flux in the total ra- 
diated power (dE/dt). When relevant, the value of ehux 
used here will be stated below. 

Figure [1] shows the dominant spectral lines from two 
relatively simple orbits, both computed to a fractional 
accuracy of Sfiux = 10^^, in the total radiated power 
{dE/dt). The orbits for these two spectra are simple 
in the sense that the motion of the test mass is very 
nearly restricted to a constant radius, and to the equa- 
torial plane of the large hole. Correspondingly, these 
spectra are also somewhat simple. For both, the peak in 
the power spectrum occurs at a frequency of 2ct;0 as one 
might guess from, for example, the waveforms computed 
by Peters and Mathews using Newtonian orbits and the 
quadrupole formula |25j. The remaining power is dis- 
tributed predominantly among three families of modes 
fixing the integer multipliers for the radial and azimuthal 
frequencies to be either zero or one. The frequencies for 
those mode families are 

{muj^ dominant power at m ~ 2 

muJ^ + ujr dominant power at m = 2 . (3.1) 
muJ^ + ujg dominant power at m = 1 

From Fig. [1] one can see that the power in any given 
mode falls off exponentially with m, at a rate that is 
determined by both the orbit geometry, and the values of 
k and n that define the family. These mode families turn 
out to dominate the spectra for all orbits with sufficiently 
small eccentricity and inclination, and the exponential 
falloff for power in modes within a fixed family turns out 
to be a general trend for these simple spectra. 

The distribution of power among modes or mode fam- 
ilies is determined by the orbital geometry. The two 
panels of Fig. [1] show that increasing orbital eccentric- 
ity draws more power into the family involving the ra- 
dial frequency, defined by modes (jl.lbp with (m, k, n) = 
(to, 0, 1). The spectrum from the system with higher or- 
bital eccentricity also has the greater number of excited 
modes which are not members of the three families that 
dominate simple orbits. The following two sections will 
demonstrate that, in general, the complexity of the spec- 
trum from any EMRI snapshot, or the number of modes 
excited by any fraction of the total power, is more sen- 
sitively dependent on eccentricity than on inclination or 
scmilatus rectum. 
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FIG. 1: The dominant spectral lines for two relatively simple orbital configurations. The parameters describing the system 
(black hole spin a, eccentricity e, semilatus rectum p, and inclination l) are shown above each panel. For the spectrum on 
the left, 80% of the total power is carried by the dominant mode and 99% is carried by the most powerful 6 modes. For the 
spectrum on the right, 77% of the total power is carried by the dominant mode, and 99% is carried by the most powerful 11 
modes. For each spectrum, 99.99% of the total power is carried by the lines shown. The red solid line traces over the tops 
of the lines carrying radiation at frequencies oimoo = rnuj^, for various values of m (which can be read off from the horizontal 
axis). Similarly, the green dash-dot line highlights frequencies ujmio = muj^ + ue (peaked at m = 1), and the dashed blue line 
shows uJmOi = muj^ + ujr (peaked at m = 2). The thicker spectral lines are members of these three mode families, while the 
thinner lines are not. 



This general rule that eccentricity governs spectral 
complexity is in accordance with the preliminary inves- 
tigation of spectral dependence on orbit geometry given 
in Ref. [l^]. There, significant waveform "voices" were 
defined by sets of frequencies ujmkn defined as follows 

azimuthal voice: fc = and rt = , (3.2a) 

polar voice: fc 7^ and 71 = , (3.2b) 

radial voice: fc = and 71 ^ , (3.2c) 

mixed voice: fc 7^ and n . (3. 2d) 

For the spectra computed in Ref. [l^l, the distribution 
of power among these voices was more strongly depen- 
dent on eccentricity than on the other orbital parame- 
ters. From the spectra in Fig. [H one might guess that 
these sets of frequencies are the best spectral classifica- 
tion scheme. The most significant of the mode families 
dominating simple spectra is exactly the azimuthal voice, 
and the other two dominant mode families are given by 
one member of either the radial or polar voices. For less 
simple orbits though, the voices defined above will prove 
a poor means of classifying spectra. For most generic 
orbit geometries, the bulk of the power is carried by the 
mixed voice, and there will prove to be a simple way of 
grouping the different members of that very large collec- 
tion of modes. 

The dominant lines of a third sample spectrum (also 
computed to a fractional accuracy of eaux = 10^^, in 
the total radiated power {dE/dt)) is shown in Fig [21 The 
orbit for this spectrum is both highly eccentric and highly 



a = 0.8M, e = 0.7, p = 6.5, i = 80 




FIG. 2: The dominant portion of the spectrum from a highly 
eccentric, and highly inclined, orbital configuration. The pa- 
rameters describing the system (black hole spin a, eccentric- 
ity e, semilatus rectum p, and inclination l) are again shown 
at the top of the plot. For this spectrum, 4% of the total 
power is carried by the dominant mode, 99% is carried by 
the top 726 modes, and about 50% is carried by the lines 
shown. Here the highlighted mode families have frequencies 
'i'lin = i^4> + i^g -\- riLOr (solid magenta line), a;o2n = 2L0g -\- nujr 
(dashed cyan line), uJ20n = 2tj0 -I- nur (dash-dot mustard- 
colored line), for various values of n. All but two of the shown 
lines are contained in these families (the two excluded lines 
are drawn slightly thinner than the others). 
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inclined. The frequencies of the dominant modes are also 
not given by Eq. (|3.ip . but are instead 

{ujff, + uje + riLOr dominant power at n = 12 
2ujs + nojr dominant power at n = 1 1 . 

2cij0 + nujr dominant power at n = 14 

(3.3) 

These mode families are not easily classified by the voices 
(|3.2p of Ref . [l^l . The values nmax of n for the dominant 
members of these mode families can be crudely approxi- 
mated (to within about 20% to 40% for the three families 
highlighted in Fig. [2|) using the conjecture (|1.6p . Equa- 
tion (|1.6|) is a good approximation to the peaks in the 
spectra derived by Peters and Mathews [2^ 



dt 



oc —< [Jfi~2{ne) - 2eJn-i(ne) 
2 

+ -Jh{ne) + 2eJn+i{ne) - Jh+2{ne)\ 
+ e^) [Jfi_2(ne) - 2Jn{ne) 



Jfl 



- — [Jn{ne)f 



(3.4) 



where here n is the multiplier of the single frequency of a 
Newtonian orbit with eccentricity e, and Jn{x) are Bessel 
functions of the first kind. Though a very crude estimator 
for the values of n describing the dominant members of 
various mode families, this formula is better than any 
other estimator that has been tried in the present work. 

By comparison to the simple spectra, the number of 
modes needed to capture any fraction of the total power 
from the complicated spectrum is much larger. This is 
demonstrated more clearly by Fig. [31 There one can 
see that, when the modes of a given spectra are sorted 
by mode number A in order of decreasing power, mode 
power falls off roughly as a power law A" , where a is orbit 
dependent, for the simple orbits. For the two simple or- 
bits, the rate of this power-law falloff is faster for the less 
eccentric orbit. If there is a power-law falloff for the more 
complicated orbit, it is not evident in the first thousand 
modes. However, given that the essentially blind algo- 
rithm from Ref. [2J| computed w 30,000 modes before 
finding the dominant thousand or so shown for the com- 
plicated orbit in Fig. [31 most of the total power is still 
captured by a surprisingly small number of modes. 

Note also that the mode which dominates the sim- 
ple orbits, with frequency 2uj^, hardly contributes to the 
spectrum of the complicated orbits. In Fig. [31 the loca- 
tion of that mode on the curve for the complicated orbit 
is indicated with a crosshairs at about A = 600, and its 
power relative to the total power « 6 x 10"^, is effectively 
insignificant. This is an extreme example of an effect that 
has been noticed in simulations of black hole binary in- 
spirals with mass ratios near unity. In the language of 
post-Newtonian descriptions of such systems (for a re- 
cent overview see Ref. [4J]) all modes other than the one 
with frequency 20;^ are "higher harmonics." The excita- 
tion of higher harmonics in those systems has been found 



iij 

V 



iiJ 

V 




FIG. 3: The normalized mode amplitudes for the orbital con- 
figurations whose corresponding spectra are shown in Figs. [T] 
and (2] sorted in order of decreasing power. The parameters 
describing the orbital geometry (eccentricity e, semilatus rec- 
tum p, and inclination t) are shown next to each curve. The 
crosshairs shown on the solid curve indicate the location of 
the mode with frequency uj = 2uj^. That same mode is the 
dominant mode for both of the other systems, as can be read 
off of Fig. (Tj The inset plot shows the same curves, but on a 
log-linear scale over a shorter range. 



to be significant in both analyt ic p arameter estimation 
studies for LISA [H, liEJilJp, IH] and fully relativistic 
numerical simulations j50l.l5lj. LISA parameter estima- 
tion studies have to date included either spin precession 
effects m, m, m Hi or higher harmonics [11, [13, E^, 
but not yet both. 



IV. SURVEY OF MANY SPECTRA 

In this section, I simulate EMRI snapshots for a large 
grid of orbital parameters and discuss how the spectral 
trends identified in the previous section vary over the 
grid. For each snapshot in the survey, the spin of the 
large black hole is taken to be a = 0.8M, and the spectra 
are computed with the code described in Ref. [IJ] using 
a requested fractional accuracy of Sfiux = 1%, m the total 
radiated power (dE/dt). The grid of orbital parameters 
is uniformly spaced in eccentricity e, inclination* l, and 
in the ratio of the semilatus rectum p to its value for 



^ It might be more natural to use a uniform distribution in cos t. 
However, this would only be a small effect on the actual values 
of inclination used. For example, the prograde orbits had ten 
inclinations uniformly distributed in t. If ten equally spaced 
values of cos t were used instead, the average difference in (, for 
any of the orbits would have been about 3° , with the maximum 
difference being about 6° . These changes would not significantly 
affect any of the conclusions in this work. 
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0.1 


3.4880 


20° 


0.1 


10.118 


110° 


0.1444 


3.7632 


25.556° 


0.1444 


10.917 


115.56° 


0.1889 


4.0388 


31.111° 


0.1889 


11.716 


121.11° 


0.2333 


4.3140 


36.667° 


0.2333 


12.514 


126.67° 


0.2778 


4.5893 


42.222° 


0.2778 


13.313 


132.22° 


0.3222 


4.8648 


47.778° 


0.3222 


14.112 


137.78° 


0.3667 


5.1401 


53.333° 


0.3667 


14.911 


143.33° 


0.4111 


5.4157 


58.889° 


0.4111 


15.710 


148.89° 


0.4556 


5.6909 


64.444° 


0.4556 


16.509 


154.44° 


0.5 


5.9662 


70° 


0.5 


17.307 


160° 


0.54 


6.2417 




0.54 


18.106 




0.58 


6.5170 




0.58 


18.905 




0.62 


6.7922 




0.62 


19.703 




0.66 


7.0678 




0.66 


20.503 




0.7 


7.3431 




0.7 


21.301 




0.74 


7.6186 




0.74 


22.101 




0.78 


7.8939 




0.78 


22.899 




0.82 


8.1691 




0.82 


23.698 




0.86 


8.4447 




0.86 


24.497 




0.9 


8.7199 




0.9 


25.295 





TABLE I: The parameters characterizing the 8000 orbit ge- 
ometries considered (all about a black hole with spin a = 
0.8M). The prograde (t < 90°) are characterized by the 4000 
possible combinations of the first three columns, and the ret- 
rograde (i. > 90°) are characterized by tlie 4000 possible com- 
binations of the Icist three columns. The two different ranges 
for p correspond to one uniform range for p/pisco- See Fig. |4] 
for a grapliical representation of the stable orbits. Note that, 
for prograde (t < 90°) orbits, pisco ~ 2.9066. Retrograde 
(t > 90°) orbits have a less relativistic pisco ~ 8.4318. 



the innermost stable circular orbit (ISCO) pisco- The 
specific values for these parameters are given in Table ID 
Of the 8000 orbit geometries shown in Table HI 728 are 
unstable. For the unstable orbits, the derivative of the 
radial potential is negative at the prescribed minimum 
radius 



d f dr 



< 



(4.1) 



where A is Mino's time parameter, related to proper time 
T hy dr = (r^ +0^ cos^ 9)dX, and where the radial poten- 
tial is 



2 2 
r 



+ {L- aEf + Q] 



Here w"^ 



,,2 



and A 



2Mr 



(4.2) 
Snapshot 



spectra were computed for each of the remaining 7272 
stable orbital configurations, represented graphically in 
Fig. m The total computational cost of simulating these 
spectra was about 2.7 CPU- years on a machine based on 
a 3.2 GHz Intel Pentium 4 Xeon processor. 
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FIG. 4: A graphical representation of the stable orbits de- 
scribed by the parameters shown in Table U] The inclination 
and minimum radius (r,nin) for the 7,272 stable orbits are 
shown. 



All of the spectra from this grid are dominated by 
mode families characterized by frequencies 



±2uj^ + 7UUr 
±U}^ + UJg + nOJr 



(4.3) 



where ± was 1 for prograde orbits and —1 for retro- 
grade orbits. To within an error on the order of 10%, 
the most dominant mode for either family had n ~ rimax, 
where rimax is given by the Peters-Mathews approxima- 
tion (jl.6p . The success of that approximation remains 
similar to that for the more complicated sample spec- 
trum in Fig. [21 and is plotted for the entire set of orbits in 
Fig. [5l Of the two dominant mode families (|4.3[) the first 
one, with {m,k,n) = (2,0,n), is the most common. The 
exceptions, dominated by the (1, 1, 7i)-family, are the or- 
bits with L nearest to 90°. This trend is demonstrated in 
Table nil The prograde orbits tend to be less easily domi- 
nated by the (1, 1, n) family. Since the orbit grid is evenly 
spaced in p/piscOi a-nd since pisco is much smaller for 
prograde orbits than for retrograde orbits, this suggests 
that the closer the orbit comes to the horizon, the harder 
it becomes for the system to channel radiation away from 
the (2, 0, n)-family and into the (1, 1, n) family. 

As with the sample orbits, the number of modes needed 
to capture the bulk of the radiated power remains reson- 
ably small. This is shown in Fig. [SI by a histogram of the 
number of modes iVgg% carrying 99% of the power. An 
example of the dependence of Ngg% on orbit geometry is 
plotted in Fig. [71 That plot shows the values of A^^gg^ 
for all the orbits with the smallest value of p/piscOi such 
that over the grid's range for the other orbital parame- 
ters, no orbits were unstable. 

It is important to emphasize that the values of Ngg% 
given in this and the following section are accurate only 
to 10%. This is because A'^ggi^ is found only after the 
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FIG. 5: The circles indicate the value of n in tJmfcn for the 
most powerful modes in the spectra from the stable orbits 
in Table HI The solid curve is the approximation (|1.6p that 
gives the peaks in the approximate spectr a (13.411 derived for 
Newtonian orbits by Peters and Mathews [25 



FIG. 6: A histogram showing the fraction of orbits, out of 
the 7,272 stable orbits in Table HI for which modes are 

required to capture 99% of the radiated power. The inset dis- 
plays the same data as a cumulative distribution (its vertical 
axis showing the fraction of orbits with Nc)q% less than the 
value on the horizontal axis). 
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20° 





25.556° 





31.111° 





36.667° 





42.222° 





47.778° 





53.333° 





58.889° 


26% 


64.444° 


94% 


70° 


100% 


110° 


100% 


115.56° 


100% 


121.11° 


100% 


126.67° 


100% 


132.22° 


1.8% 


137.78° 


0.3% 


143.33° 


0.3% 


148.89° 


0.3% 


154.44° 


0.3% 


160° 


0.3% 



TABLE II: The fraction of orbits, from the grid of orbits 
described in Table U dominated by a mode with frequency 
^mkn, where {m,k,n) — (±1,1, n), as a function of incli- 
nation L. The remaining orbits are dominated by modes 
with (ra,k,n) = (±2,0, n). Here ± is 1 for prograde orbits 
(t < 90°) and -1 for retrograde orbits (t > 90°). Note that 
since the grid of orbits is uniformly spaced in p/pisco, rather 
than in p, the prograde orbits have smaller values for p than 
do the retrograde orbits. 
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FIG. 7: The number of modes required to capture 99% of 
the radiated power, A^g9%, as a function of eccentricity e and 
inclination t for a fixed value of p/pisco- Since pisco is larger 
for retrograde orbits, the top block of this plot corresponds 
to orbits with larger values of p. 



algorithm from Ref. [24| computes a much larger number 
of modes N in an effort to determine the total power 
to its requested fractional accuracy of Eaux ~ 1%- Once 
that algorithm terminates, the modes that it computed 
are sorted in order of decreasing power such that 



dE\ _ A /dE\ 



(4.4) 



where (dE/dt)^ decreases with A. The value of iVg9% is 
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the smallest possible value satisfying 



10" 



,x 10 



TV 

E 

A=Ngg' 



dE\ 



<,0.„l,(f 



(4,5) 



Spot checking its dependence on efiux for a few sample 
spectra gave an expected accuracy ~ 10%. 

Figure [7] also shows that, as with the sample orbits, 
iVgg% is more strongly dependent on eccentricity than on 
inclination and semilatus rectum. This is especially sig- 
nificant for the prospect of observing intermediate mass 
ratio inspirals with ground-based detectors like LIGO, 
since the systems that have been estimated to be the 
most likely candidates for being observed are those with 
especially small eccentricity, typically e < 10^^ and at 
most e w 0.1 [l3l. For LIGO, it is likely that waveforms 
needed for detection need only contain a few to ~ 10 
modes. This both simplifies LIGO's task of detection, 
since the waveform snapshots will not be very compli- 
cated, and hardens its task of spacetime mapping, since 
correspondingly less information will be observable. 



V. SPECTRA FROM A KLUDGED INSPIRAL 

In this section, I describe how general spectral charac- 
teristics should be expected to evolve during an EMRI. 
The snapshots studied here will be sampled at approx- 
imately 12-hour intervals over the final three years of a 
single kludged EMRI thought to be typical of the kind 
that could be observed by LISA. The kludged trajec- 
tory through orbital parameter space e(t), L{t), and pit), 
was provided by Jonathan Gair, and was numerically 
computed according to the prescription introduced in 
Ref. [56| . Their method for approximating the trajec- 
tory is based on an eclectic combination of approxima- 
tions including post-Newtonian equations for the radia- 
tive fluxes of energy and angular momentum, numerical 
fits to Teukolsky-based calculations of fluxes for circular 
and equatorial orbits, as well as some uncontrolled ap- 
proximations for the evolution of Carter's constant. Since 
only power spectra will be discussed here, the results are 
independent of any evolution for the positional orbital 
elements, or Xmkn hi Eq. (|2.28p . While one would expect 
minimal accuracy from such an array of approximations, 
these kludged trajectories have been shown to exhibit 
a stunning degree of agreement with more accurate cal- 
culations. The integrated overlap between approximate 
waveforms based on the kludged orbit trajectories and 
waveforms constructed from black hole perturbation the- 
ory alone is often about 95% [13]. So a kludged orbit 
trajectory is likely more than adequate for the present 
purpose, since the snapshots themselves will still be ac- 
curate up to leading order in the mass ratio, and since 
for the majority of the spectra examined in the previ- 
ous section, the general character of a spectrum does not 
change dramatically for small perturbations to the orbital 
parameters. 





1.5 2 
t (years) 



FIG. 8: The kludged evolution for a selection of the orbits 
principal elements during an inspiral of the sort that could be 
observed with LISA. For this system, the large black hole has 
mass M — 10^ A/©, and the magnitude of its spin angular mo- 
mentum is a = 0.9A/. The smaller black hole is nonspinning 
and has mass ^ = lOM©. The initial orbit has eccentricity 
e = 0.79, semilatus rectum p = 8.8, and inclination l — 43°. 
The final orbit has eccentricity e — 0.10, semilatus rectum 
p = 3.2, and inclination t — 45°. The orbit trajectory e{t), 
p{t), and L{t) for this inspiral was computed by Jonathan 
Gair using the algorithm described in Ref. [s^. The approx- 
imations used to compute the trajectory should become less 
accurate with time as the system becomes increasingly rel- 
ativistic. The top panel shows the evolution of the orbital 
frequencies fr,e,4, ~ ujr,e,<t,/ {'2t^) with an inset of the same but 
over a different range. The bottom panel shows the evolution 
of orbital eccentricity e, with an inset of the same but over a 
different range. 



The kludged orbit trajectory for the inspiral examined 
here is shown in Fig. [8p The principal orbital elements 



There is no special reason for using a black hole spin of a = 0.9M 
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FIG. 9: The initial (blue, or dark grey) and final (red, or 
light grey) power spectra for the inspiral shown in Fig. [8] An 
animated movie of the complete spectral evolution is available 
upon request. 
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FIG. 10: The evolution of the number of modes needed to 
capture 99% of the radiated power Ngg% for the inspiral shown 
in Fig. |5] The inset shows a close-up view over about the last 
week of the inspiral. 



evolve slowly and smoothly throughout the majority of 
the InspiraL In the final few days of the inspiral, where 
the trajectory should be least accurate and where the 
adiabatic approximation itself should begin to fail, the 
eccentricity begins to rise with time, and the orbital fre- 
quencies rapidly diverge from each other. The trajectory 
ends when it has evolved onto an unstable orbit, at which 
point the binary would merge to form a single black hole. 

Snapshot spectra were computed at roughly 12-hour 
intervals along the orbit trajectory in Fig. [51 As was done 
for the grid of spectra in the previous section, these spec- 
tra were computed with the code described in Ref. [l^l 
using a requested fractional accuracy of esux = 1%, hi 
the total power (dE/dt). There were 2,143 snapshots in 
all, and the total computational cost of simulating them 
was about 1.5 CPU-years on a machine based on a 3.2 
GHz Intel Pentium 4 Xeon processor. The initial and 
final snapshot spectra from this sequence are shown in 
Fig. [5) The inspiral's initial spectrum is similar in char- 
acter to the one shown in Fig. [21 due to the large initial 
orbital eccentricity. The final spectrum is more similar 
to the ones in Fig. [TJ The dominant mode families for 
both the initial and final spectra have (m, fc, n) = (2, 0, n) 
and (1, 1, n). Those two families make up the two largest 
arcs or lobes of lines in the initial spectrum. In the fi- 
nal spectrum, the lobes are much more narrow, and are 
not as easily identified by eye. For the final spectrum, a 
more efficient definition of mode families might instead 
use m as the free index, rather than n. For example, 



in the kludged trajectory as opposed to a = 0.8M, as is used for 
the grid of orbits in the previous section. The difference is due 
to the trajectory and the present work not being produced in 
parallel. 



the strongest lines along the upper right edge of the final 
spectrum are (m, 0, to), peaked at to = 2. 

Figure [10] shows the evolution in spectral complexity 
during the inspiral. As the binary becomes more cir- 
cular, the number of lines carrying 99% of the spectral 
power A^ggi^ decreases with time until the last few days, 
at which point both the kludge and adiabatic approx- 
imations should fail. Over those last few days the or- 
bital eccentricity rises and Ngg% jumps from about 125 
to 250. This effect is typical of both kludged trajecto- 
ries and the more accurate Teukolsky-based trajectories. 
While it may be a physical effect, it always occurs sus- 
piciously in regimes where the adiabatic approximation 
and kludge should be least accurate. The question of 
whether or not it is physical might best be addressed by 
the numerical relativity community. The causes of the 
smaller abrupt jumps in Ngg% (for example, at either 
end of a six month gap centered at about one year) is 
unknown, however, they are within the expected accu- 
racy ~ 10% whereas the final jump from 125 to 250 is 
not. Repeating the snapshot calculations with a smaller 
requested overall accuracy in the total power eaux would 
likely eliminate most, if not all, of these small jumps. 



VI. VERIFICATION OF BLACK HOLE ORBITS 

In this section I suggest EMRI detection algorithms 
based on the spectral trends found above and discuss the 
scientific meaning of successful detections. The tone of 
the discussion is meant to be exploratory rather than 
exhaustive. That is, it is meant to outline the sorts of 
data analysis strategies that the simulations described 
above suggest would be useful, rather than to look in 
great detail at any one algorithm. 

This section will be divided up into two subsections. 
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The first subsection outlines algorithms that could search 
for systems without radiation reaction. These are really 
just of interest as building blocks for more complicated 
searches since only a small subset of EMRIs would be 
observable without accounting for the influence of radi- 
ation on the binary. The second subsection describes 
how these building blocks would be used to search for 
realistic systems affected by radiation reaction and gives 
a detector-independent estimate of how well these algo- 
rithms might perform in the best of circumstances. 



A. Without radiation reaction 

For times that are sufficiently short for the orbit of the 
captured mass to be unaffected by radiation^, the gen- 
eral expression for EMRI snapshots (jl.ip is accurate up 
to corrections that are at most of order jj? /M'^ and are 
perhaps even as small as ^ /i^/M^ with the principal 
frequencies adjusted according to the second order met- 
ric perturbation [43j . In this subsection, I will outline a 
detection strategy for a signal of this form. Truncating 
the general expression for an EMRI-snapshot waveform 
(jl.ip to a finite number of N modes, rewriting it as a 
single sum, and explicitly showing the positional phase 
elements, gives the following 



N 



exp{-i[mAW0(i - t^) 

A=l 

-I- kMJe{t -te) + n^ujrit - tr)]] ■ (6.1) 

I wish to consider the prospects of using this truncated 
expression (j6.ip explicitly as a phenomenological tem- 
plate. 

By phenomenological templates, I mean templates 
for which some of the free and measurable parame- 
ters will have no immediately obvious physical meaning. 
Traditional, or nonphenomenological templates are con- 
structed as follows. You start by declaring the 17 param- 
eters which completely determine the waveform (e.g. the 
two masses, and 3-vectors for the position and velocity 
of each object as well as for the spin of the larger ob- 
ject). You then solve Einstein's equation or some per- 
turbed form of it to completely determine ft.+.x- This 
last step is equivalent to determining all of the frequen- 
cies Wr,fl,0, phase shifts tr_e,tj>, and Fourier coefficients /ia 
in the model (|l.ip . The phenomenological templates are 
computed in a much simpler way. You first chose values 
for the frequencies, phase shifts, and Fourier coefficients, 
and then you simply evaluate the sum over modes (|6.1|) . 
For these waveforms, the frequencies, phase shifts, and 



Fourier coefficients are themselves the free quantities to 
be measured. For a phenomenological template with N 
modes, there are then 5A^ -I- 6 free parameters. Of those 
there are 2iV-|-6 real numbers (the three frequencies, the 
three phase shifts, and the N complex mode amplitudes), 
as well as 3N integers (the frequency multipliers m\, k\, 
and tia). 

The most obvious objection to the use of these phe- 
nomenological templates in gravitational wave searches 
is the large dimensionality of the parameter space. For 
templates that allow for any more than N = 2 modes, 
you will have more free parameters than the traditional 
templates^, and of those only the three frequencies uj^.e^r 
will carry immediate physical meaning. In the traditional 
scheme, however, the step of going from the complete set 
of 17 parameters to the waveform is extremely costly. 
Since both template families are equally good matches 
to the true waveforms, and since either way you will be 
dealing with a significant number of template parame- 
ters, it may be worth adding dimensionality to the tem- 
plate parameter space in exchange for not having to solve 
Einstein's equation. 

Though the large dimensionality may seem daunt- 
ing, especially for those familiar with efforts focused on 
sources with just a few degrees of freedom, it is not unre- 
alistic for elaborate algorithms to identify signals charac- 
terized by a large number of parameters. For algorithms 
that match against templates over a sequence of increas- 
ingly dense grids on the model parameter space, com- 
putational cost grows as a power law where the power 
is proportional to the number of free parameters. Algo- 
rithms designed for more complex models have costs that 
instead grow only linearly with the number of parame- 
ters. Algorithms of this nature have been used in the 
mock LISA data challenges [l4| to recover ~ 10"* white 
dwarf binaries. 

I now discuss how well phenomenological detection al- 
gorithms can be expected to perform in the best case sce- 
nario where the algorithm has no difficulty in selecting 
the correct parameter values. This can be done without 
reference to specific instruments, and in a sky averaged 
sense, by studying the distribution of the power among 
various modes relative to the total power radiated by the 
binary. More sophisticated estimates that account for de- 
tector characteristics and variation of signal parameters 
arc certainly possible, but they are beyond the scope of 
this paper. 

For the task of searching only for EMRI snapshots with 
phenomenological templates, the example spectra from 



For a subset of observable systems |22 | these "short" times are 
actually longer than the longest amounts of time for which fully 
coherent search algorithms arc computationally affordable (a few 
weeks @|). 



^ The specific value N = 2 here is not exact. The three integer 
parameters are much simpler degrees of freedom from the stand- 
point of a search (discrete and reasonably confined) . In the same 
vein though, not all of the 17 traditional parameters are intrin- 
sic. So while the value of N for which both the traditional and 
phenomenological template spaces are dimensionally equivalent 
is surely small, the exact value is not obvious. 
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Sec. mil suggest estimates of the minimum scale of the 
phase space dimensionahty. For example, systems with 
eccentricities of about 1% will produce spectra similar 
to Fig. [TJ So a ground-based phenomenological search 
for these systems would require < 10 modes, result- 
ing in a model with 56 parameters (26 real numbers and 
30 relatively small integers). For snapshots thought to 
be typical of EMRIs that could be observed with space- 
based detectors, one would need N ~ 10^ to lO'^. This 
would mean a model with about 500 to 5000 free param- 
eters, still far fewer than what has been demonstrated 
already for white dwarf binaries [lij . 

If a search for phenomenological EMRI snapshots were 
successful, only the general waveform model would be 
verified. One would only be able to claim detection of 
some signal with a discrete triperiodic spectrum with 
some measured fundamental frequencies uj^^g^r^ since the 
physical meaning of the other unknown parameters is 
convoluted. In this event, one could then turn to the 
underlying physical model. Finding a set of its 17 free 
parameters that best reproduce the parameters measured 
in the phenomenological search would then confirm the 
physical model to some more explicit level of uncertainty. 
Failing to find parameters for the underlying physical 
model might mean that the snapshot was produced by 
a test mass moving along a geodesic of some non-Kerr 
spacetime, since many (but not all) candidates for such 
orbits are also triperiodic [H, [5§| . 

Some intermediate level of model verification is also 
possible and could reduce the dimensionality of the pa- 
rameter space for the phenomenological templates. For 
example, sixteen mode families are needed to capture 
99% of the power radiated by the initial snapshot from 
the inspiral examined in the previous section. The distri- 
bution of power among these families is shown in Table 
mil As is true for all the snapshots simulated in this 
paper, this one is dominated by the mode families with 
(m, k, n) = (2, 0, n) and (1, 1, n). And as is typical, those 
two mode families carry most of the power, 73% here. 
In an effort to simplify the phenomenological waveform 
model, one might restrict it to include only those mode 
families. This specific model would eliminate (ttja, kpC) 
from the template parameter space by fixing them to ci- 
ther (2, 0) or (1, 1), and would create two new parameters 
specifying the number of modes in each of the two fam- 
ilies. This would reduce the number of free parameters 
from 57V 6 to 37V -I- 8. 



B. With radiation reaction 

The scenario describe in the previous subsection is only 
immediately useful for EMRI's with the most extreme 
mass ratios. There is no compelling reason to expect 
those systems to be especially common, or to even con- 
sider them reasonable targets at all. The purpose of 
studying these simple systems is to construct from the 
results an approximate description of more generic EM- 
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TABLE III: The 16 mode families needed to capture 99% of 
the power radiated during tlie first 12 hours of the inspiral 
discussed in Sec. |V] Each family is defined by fixing m and 
k. The most powerful member of each family has frequency 
muj^ + kuje + Wmaxt^'r ■ The ratio of the power radiated by each 
family to the total power from this first 12 hours is given by 

(E) /(E). For this system, the spin of the black hole 

\ / mk \ ' 

is a = 0.9A'f, and the orbit of the test mass has eccentricity 
e = 0.79, semilatus rectum p = 8.8, and incUnation t = 43°. 



RIs that respond to their own radiation. That is, we wish 
to describe the radiation of a generic adiabatic EMRI as 
a slowly evolving sequence of EMRI snapshots. Here I 
will now outline how the phenomenological templates for 
EMRI snapshots could be modified to describe the more 
general class of adiabatic EMRIs. There arc many ways 
that this could be done. Although a detailed study of 
specific models would be valuable, it is beyond the scope 
of this paper. Instead, I aim to be as general as possible 
and will steer away from discussing any specific imple- 
mentation. 

For an adiabatic EMRI, the motion of the small object 
is described by a solution of the geodesic equation for the 
Kerr spacetime, but with the orbital elements replaced 
by quantitates that evolve slowly. The Teukolsky equa- 
tion can provide the leading order radiative changes to 
those quantities, and more sophisticated techniques are 
envisioned for describing both conservative and radiative 
effects. In the spirit of trading calculation difficulty for 
added dimensionality, every quantity that was constant 
for the snapshot model (|6.ip could in principle be re- 
placed by simple, one or two parameter models. 

To illustrate this, consider the orbital frequencies 
for i = r, 0, (/). These can be taken to drift linearly with 
time 

LOi — > uji + ujit . (6.2) 
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where I have introduced new constants uji oc fj,/M. This 
simple model fits the first year of the three frequency 
trajectories shown in Fig. [8] with an average fractional 
accuracy of about 1%. Other models can of course do 
better. For example, Peters and Mathews derived an 
expression for radius as a function of time in the case of 
slow circular inspirals. Combining that result, Eq. (5.9) 
of Ref. nil, with Kepler's law Muj = (r/M)-^/^, gives 



Mlu 



MI 



256 fi t 
~5~MM 



-3/8 



(6.3) 



For the more general case of fast generic motion, one 
might want to try a model with a similar form 



(6.4) 



This model fits the first year of the three frequency tra- 
jectories shown in Fig. [8] with 7, w 1.6 and an average 
fractional accuracy on the order of 10~^. It performs 
similarly at later times, but not with the same values of 
the parameters a,, /3j, and 7^. Unlike the simple linear 
model however, this one has three parameters instead of 
two. 

To complete the construction of phcnomcnological 
templates for adiabatic EMRIs, similar models can be 
concocted for the other parameters of the snapshot tem- 
plates (j6.ip . the complex mode amplitudes, and the posi- 
tional phase elements. If one-parameter models like the 
linear frequency drift are used, then the new template is 
given by Eq. (|6.ip with 
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again for i = r, 9, (f>. For this case, the dimensionality of 
the template parameter space is doubled to lOA^ -I- 12. 
Of those, only the 3N frequency multipliers are inte- 
gers. These phenomenological templates are similar in 
nature to the time-frequency search methods which have 
been successfully demonstrated in the mock LISA data 
challenge [63|. They differ from those in that they can 
accommodate coherent integration and could also more 
naturally include specific schemes for evolving both the 
principle and positional orbital elements. 

Given that EMRI snapshots tend only to be dominated 
by a small number of mode families, it is likely that a 
slightly simpler waveform model could be used for adi- 
abatic EMRIs. Figure [11] demonstrates how successful 
two such hypothetical detection algorithms might be if 
searching for the kludged inspiral from Sec. El If one 
were to search for this waveform using a model that in- 
cluded only the dominant three mode families from the 
sample complicated snapshot in Fig. [2] one would recover 
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FIG. 11: A time-dependent measure of success for two hypo- 
thetical detection algorithms searching for the inspiral from 
Fig. [S] The algorithm for the dashed blue curve captures the 
power radiated at all frequencies of the form uj^kn, where 
{m,k,n) = (1,1, n), (0,2, 71), and (2,0,n), for any n. The 
solid curve captures power radiated at all frequencies in the 
16 mode families needed to capture 99% of the initial power, 
where here a mode family is defined as all modes with fre- 
quencies ijJrnkn for some fixed m and k. 



71% of the inspiral's total power. Assuming this could be 
done with one-parameter models for the adiabatic evolu- 
tion of the waveforms parameters, such a model would 
have 6A^ -I- 18 free parameters (by eliminating 2N fre- 
quency multipliers and their 2N linear drifts, and by 
adding three new integers specifying how many modes 
are in each family, as well as their three linear drifts). 
For this example EMRI, this scheme included N = 121 
modes. So 714 free parameters would have been needed 
to recover 71% of the power. A similar hypothetical al- 
gorithm which captured all the power carried by the 16 
mode families that make up most of that inspiral's ini- 
tial spectrum (Table IIII[) would recover 98% of the total 
power with 6A'^ -|- 44 free parameters (by eliminating 2N 
frequency multipliers and their 2A^ linear drifts, and by 
adding 16 new integers specifying how many modes are 
in each family, as well as their 16 linear drifts). For this 
example EMRI, this scheme included = 411 modes. So 
2,510 free parameters would have been needed to recover 
98% of the power. 

It should be emphasized that any gravitational wave 
detections following from the use of these phenomenolog- 
ical waveform models would not necessarily yield either 
the parameters that completely determine an EMRI (po- 
sition, masses, etc) or the spacetime map that general rel- 
ativity predicts is encoded in the radiation. They would 
however verify the detection of waveforms predicted to be 
produced when a test mass perturbs a rotating black hole 
by moving through an adiabatic sequence of its bound 
geodesic orbits. They would also measure the evolution 
of the three fundamental frequencies throughout that se- 
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quence, or equivalently the evolution for any other set 
of principle orbital elements. It is possible that restrict- 
ing phenomenological EMRI waveforms to include only 
the mode families that are most commonly dominant in 
the snapshots simulated here may alone be enough of 
a constraint to keep EMRIs into non-Kerr black hole 
candidates from triggering a detection. However, with- 
out more work with the snapshots from such EMRIs one 
could not say so with any certainty. One would have to 
be content only to have verified that radiation from an 
adiabatic sequence of black hole orbits could have trig- 
gered the detection. 



VII. CONCLUSION 

The number of significant modes in generic EMRI- 
snapshot spectra has been shown generally to be much 
more manageable than one might have guessed from ear- 
lier truncation algorithms [l^]. This should lead to im- 
proved truncation algorithms, which will reduce the cost 
of future data analysis efforts. Such improvements should 
exploit the trends observed here in the relationship be- 
tween orbit geometry and spectral signature. The ability 
to predict the multiplier n « Timax of the radial frequency 
for the dominant modes using a formula based on such 
simple approximations [1^ is encouraging. It suggests 
that many of the trends in these spectra might be under- 
stood analytically using more recent tools [6ll, [gI, [6^ . 

The detection algorithms that are suggested here for 
verifying minimal aspects of relativity and black hole 
physics may ultimately be used in future gravitational 
wave detections. However, more work is needed to deter- 
mine whether or not they are cost-efficient and science- 
efficient alternatives to traditional search techniques. 
Another possibly interesting area for future work is to 
explore the dependence of the snapshot spectra on the 
spin of the larger black hole. It has been implicitly as- 
sumed that the large values of spin considered are some- 
how representative of observable EMRIs. This is reason- 
able since the few existing measurements due to modeling 
x-ray spectra from galactic nuclei suggest near-maximal 
spins [U, [1^ . Still, future work that tests the general- 
ity of these parameter values by simulating other EMRIs 
would be worthwhile. 
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APPENDIX A: BI PERIODIC FORM OF 
TIME-DEPENDENT SPATIAL COORDINATES 



Here I derive the bi-periodic form (|2.3p of the spa- 
tial Boyer-Lindquest coordinates r, 9, and (jj for bound 
geodesies function of the time coordinate t. 

The bi-periodic forms of both the radial and polar co- 
ordinates follows immediately from Sec. IV of Ref. 2l| by 
replacing f[r{t), 9{t)] with r{t) and 6{t). The resulting re- 
lationships between the coefficients in the Mino-time se- 
ries expansion (j2.ip and the coordinate-time expansions 
are given by 



(Al) 



TrTe 
(27r)2r 
dt 



dXr 



27r/Te 



dXf) 



(A2) 



where r(Ar) and 9{Xe) are given by their Mino-time scries 
expansions with A = A,, and A = Ae, respectively. 



and where 



StiXr, Ae) = y^Jkr, 



-ikT e\e—inT 



(A3) 



(A4) 



The A-derivative of t can be evaluated analytically as 
given by Carter's first order geodesic equations, or it can 
be found from 



^[r{Xr),6{Xe)]^T 



■ ^(fcTe + nTr)tkne~ 



■ikT e^e —itiTrXr 



kn 



(A5) 



The derivation of the bi periodic form for (/)(t) is 
sfightly different. First, write t{X) as 



t = VX'V5t . 



(A6) 
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Multiplying by w^, and rearranging terms gives 

T^A = n^t - n^st , (A7) 

since to^T = T^. Now write ^(A) as 

cl>^T^X + Scj) , (A8) 

and insert the above expression (|A7p for the first term to 
find 



where 



= 5(f> — ui^St . 



(A9) 



Treating (50 as a function of r and 6 now allows it to 
be used in place of /[r(A), 61(A)] in Sec. IV of [2l|. This 



gives the following expression for the coefficients of the 
bi periodic form of 4){t): 



dt 



X — [r(A^),6l(Ae)]e^('='^''+""'-^''*^^'-^^''^ 



X S(l}{Xr,Xe)e 



ikT Xe -\-inT r Xr 



(All) 



where 



(AlO) (50(Ar, Ag) = y^((/)fc71 ~ UJ^tkn)e~ 
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